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INTRODUCTION 


The  purpose  of  the  report  was  to  determine  the  effect 
of  sensor  errors  in  the  line -of- sight  rate  on  the  miss 
distance  for  an  accelerating  target.  Since  ARRADCOM  is 
responsible  for  developing  fire  control  systems  for 
antitank  application,  it  was  decided  to  develop  a 
methodology  which  would  accept  the  basic  error  per 
measurement;  and  then  by  using  a  number  of  these  measurements 
in  succession  to  develop  a  smoothing  and  prediction  scheme 
which  would  quantify  important  parameters  at  the  time  of 
launch.  These  estimated  parameters  at  the  time  of  fire 
are  important  in  determining  the  gun  lead  angle  and  the 
resulting  miss  distance  at  the  target.  Also,  measurements 
taken  by  the  sensors  are  not  perfect  and  result  in  errors 
in  the  gun  lead  angle  and  therefore  in  miss  distance. 

BACKGROUND 


Previous  studies  have  concentrated  on  various  filter 
models  and  classes  of  moving  vehicle  using  actual  target 
dynamic  data.  Continuous  tracking  was  assumed  with  no 
parametric  sensitivity  analysis  to  show  impact  errors  as  a 
function  of  short  tracking  times  and  sensor  noise  (ref  1-3) . 
This  analysis  presents  work  done  from  January  through  March 


of  1979. 


METHODOLOGY 


The  design  methodology  is  shown  in  figure  1. 


Figure  1.  Target  and  launcher  geometry  at  time  of  fire. 


Since  0  is  the  Y  component  of  the  velocity  of  the  target 
divided  by  range  X, 


0  =  VT  •  SIN  a  /X 


(1) 


It  is  assumed  that  the  target  is  following  a  path  in 
the  prelaunch  phase  of  0  which  equals  a  constant  prior  to 
firing.  As  a  result,  0  will  vary  linearly  with  time.  A 
least  squares  fit  of  the  data  is  therefore  a  best 
approximation  by  drawing  a  straight  line  through  the 
measured  data  points  in  the  prelaunch  phase.  Employing 
linear  regression  techniques,  one  obtains  (see  appendix  A) 


E0  = 


1 

Det 


NMAX 

Z  Cn  •  $Nn 
n=l _ 

3 


(2) 


Where  Cn  is  defined  by  (appendix  A) 


<NMAX  \  /NMAX\  / 

E  tn 2J~  tn  '^Z  tn^-  ^tn  •  NMAX 


n 


NMAX  \ 

E  tn  )  •  TF 
n=l  /  (3) 


and  where  the  summation  signs  are  from  nml  to  nmNMAX. 
Therefore, 


NMAX 

Det  =  NMAX  •  E  tn2  - 
n=l 

One  then  obtains  solutions  for  the  next  two  equations 
(see  appendix  A) . 

The  first  equation  states  that  at  intercept,  the  X 
distance  that  the  projectile  travels  away  from  the  launcher, 
TINT  •  VP  •  COS  AP,  must  equal  the  distance  of  the  target 
form  the  launcher.  This  distance  is  the  initial  distance 
at  time  of  fire,  Xo ,  plus  the  distance  travelled  due  to  the 
initial  velocity  in  the  X  direction  TINT  *  VT  '  COS  a ,  and 
the  distance  travelled  due  to  the  component  of  target 
acceleration,  A  COS  a'  %  TINT2  and  results  in 


NMAX 

a  tn) 

n=l 


(4) 


TINT  •  VP  •  COS  AP  =  X0  +  VT  .  COS  a  •  TINT 

+%  (A  COS  a*)  •  TINT2  (5) 

Since  at  intercept  the  distance  travelled  by  the 
projectile  in  the  y  direction  must  equal  the  distance 
travelled  by  the  target  in  the  y  direction,  one  obtains 


TINT  •  VP  •  SIN  AP  *  VT  •  SIN  a  •  TINT 

+%A  •  SIN  a’  •  TINT2  (6) 

Equations  5  and  6  are  then  solved  for  TINT  and  AP. 

One  uses  the  nominal  intercept  time,  TINT,  to  calculate 
the  miss  distance.  These  subsequent  calculations  are  based 
upon  the  launcher  measurements  and  measurement  errors . 

The  line- of- sight  rate  error  is  the  main  focus  of  interest. 
Therefore,  errors  arising  from  parameters,  other  than  the 
line  of  sight  rate  error,  are  not  considered  here.  The 
influence  of  the  error  on  transverse  velocity  and 
acceleration  estimates  are,  however,  assessed,  since  poor 
estimates  will  result  in  an  increased  miss  distance  at  the 
target.  For  example. 


HD  =  TINT  •  (VTJ  -  VTJP) 

+  <™112  .  (AJ  .  AJP) 

By  definition  of  Ee  (appendix  A) 
Ee  =  (9P  -  S)/fi 

and 

Ee  -  (eP  -  e)/e 

since 

fl  -  VTJ/X0 

and 

6P  -  VTJP/Xq 


(7) 

(8) 

(9) 

(10) 

(ID 
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substituting  (10)  and  (11)  into  (8) 


♦ 


VTJP  -  VTJ 

Ee  =  /e  (12) 

or 

Ee  =  (VTJP  -  VTJ)/(Xq  •  6)  (13) 

or 

VTJP  -  VTJ  =  Ee  '  XQ  •  0  (14) 

From  Greenwood,  page  37,  equation  2-23  (ref.  4),  ex  =  i 
and  ei  =  j.  Solving  the  j  component  for  V.  i.e.,  0  in 

the  nomenclature  of  this  report 


0  =  (AJ  -  2  •  0  •  VTI)/X0 
and  similarly 


(15) 


0 P  =  (AJP  -  20 P  •  VTI)/XQ 

Substituting  (15)  and  (16)  in  (9)  yields 

(AJP  -  2eP  •  VTI)  -  (AJ  -  20  •  VTI) 

Ee  “  Xq0 

or 

(AJP  -  AJ)  +  2 VTI  (0  -OP) 

Ee  -  ( Xd'9) 

Substituting  (8)  into  (18),  yields 

E0  =  [(AJP  -  AJ)  +  2  •  VTI(-E0  •  0)]/  (Xo  *  0) 

From  (19)  and  solving  for  AJ  -  AJP, 

AJ  -  AJP  =  -Ee  •  Xc  ■  e  +  2  •  VTK-EO  •  0) 

From  (13) , 

VTJ  -  VTJP  =  -E0  •  XQ  •  0 
Substituting  (20)  and  (21)  into  (7) ,  one  finds 


(16) 

(17) 

(18) 

(19) 

(20) 

(21) 


t 


t 


4 


TINT2 

MD  =  TINT  •  (-E0  •  XQ  •  §)  +  2 

(-E0  •  XQ  0)  +  2  •  VTI(-E§  •  S) 


or 


MD  =  E0 • 

+  E0 


TINT2 

-TINT  •  XQ  0  -  2 


2  *  VTI  •  6 


TINT2 


or 


"2  ‘  x0  0 


TINT2 


Inserting  (25)  and  (26)  into  (24)  yields 
MD  =  E0  *  Ki  +  E0  *  K2 


(22) 


(23) 


MD  = 

E5 1 (TINT  •  XQ 

+  T~ 

2  •  VTI)  (-1)5 

..  [  TINT2 

+  E0  1  -  2~ 

X 

O 

CD 

_ 

(24) 

(24) 

let 

r 

TINT2 

•  ^ 

(TINT  •  XQ  + 

“ 2~ 

•  2  •  VTI)  (-1)  0 

(25) 

"  TINT2 

- 

K2  = 

-  2  *  X0 

0 

(26) 

(27) 


Using  linear  regression  techniques  (see  appendix  B) 


NMAX 

^  .  _±_  E 

E0  =  q  1  Det  n=l  0N 


n 


n 


NMAX 
NMAX  Z  t. 


n=l 


n 


(28) 


Let 


CTDD 


.-k 

L 


NMAX 

NMAX  -  Z  t 
n  -  1 


n 


(29) 


Inserting  (29)  into  (28)  yields 


NMAX 

..  ±  1  z 

E0  “  e  •  Det  n=l 


0Nn  *  CTDDn 


(30) 


(2)  and  (30)  in  (27) 
NMAX 

z 

Ki  n=l 
HD  =  Eet 


Cn  '  8» 


n 


K2 
+  — 


6  •  Det 


NMAX 

z  e 

n=l 


N 


•  CTDD 


n 


n 


(31) 


or 


NMAX  .  /Cn  •  K 
n=l  “ 


MD  -  1  -  V  Det  •  6 


1  +  CTDDn  ■  K2  , 
e  Det  j 


NMAX 


(SIGMA  MD)2  =  Z  a2  6n. 


SIGMA  MD 


n=l 


NMAX 


,  /' cn  ,Ki 

[n  Det  •  0  0  •  Det 


+ 


(32) 


CTDDn  •  K21 


JL  JLJJLJ-q  *  is.2 

•Det  / 


(33) 


'  n=l 


.  (Cn  ■  Ki  ^  CTDDnK2\2'  %  (34) 
0 Nnl  Det  •  e  0  •  Dety 


If  0  2 ©M  is  constant  with  n,  by  substituting  °20M,  one 
obtains1* 


SIGMA  MD  =„  0 


a  °N 


N  MAX  /Cn  •  ¥LX  CTDDn  •  K2' 


2  -% 


E  V 

n-1  V 


Det  *0  0  •  Det 


(35) 


% 


# 
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CONCLUSIONS 


By  using  this  approach,  it  has  been  possible  to  present 
*  the  main  statistical  concepts  and  to  relate  the  variance  in 

the  tracking  rate  noise  to  the  variance  in  the  mil  error 
at  the  target.  This  concept  takes  into  account  the  smoothing, 
prediction,  projectile  flight,  and  a  second  order  prediction 
scheme.  Effects  of  other  important  variables  such  as  data 
rate,  smoothing  time,  firing  delay  (due  to  remaining 
computations  and  final  aiming  corrections  required  after  the 
last  measurement) ,  target  velocity,  target  velocity  direction, 
target  acceleration,  target  acceleration  direction, 
projectile  velocity,  range  and  gun  lead  angle  have  also  been 
included . 


RECOMMENDATIONS 

Analyses  should  be  performed  which  show  the  effects  of 
the  above  parameters.  Further  recommendations  will  be 
made  pending  completion  of  these  parametric  investigations. 
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APPENDIX  A.  DERIVATION  OF  ERROR  IN  LINE  - 
OF -SIGHT  RATE 


V 


» 


» 


Figure  A-l.  Smoothing  process, 
=0o+0o  •  tn  +  eNn 


(A-l) 


e 


en 


(A-2) 


(where  tn  =  Dt  •  (n-1)) 

H  =  measured  value 
e  =  Estimated  value 
N  =  noise  error' 
n  =  n-th  measurement 

For  best  fit  of  dash  line  to  the  5™  pointSj  the  sum  of  the 
differences,  squared 


nx 

=  S  -  i  (eM  -  een)2 
n=l 


(A- 3) 


By  substituting  (A-l) 


and  (A-2)  into  (A-3) ,  one  has 


P 


* 


8 


nx 

S  =  E 
n=l 


+  <«<,  -  Oo>cn  -  «Nn 


(A- 4) 


To  minimize  S  with  respect  to  the  two  variables,  0Qe  and 
e  oe>  cake 


„  nx 

££  =  o  =  z 


ae 


oe 


n^l 


+  <«oe  -  «0)<4>  -  8Nn 


(A-5) 


3S 


nx 


££  -  0  =  z  2 


30. 


n=l 


<0O-  -  0o) 


+  <eoe  "  e0)t:n  -  ®Nn 


~n 


or 


(A-6) 


nx 

z 

n=l 


<5oe  -  8o>  +  (9oe  -  ycn 


nx 

=  Z  .  ®Nn 
n=l 


(A-7) 


nx 

z 

n=l 


(0  ^  -  0  )  +  (0rt  -  0  )  t-. 
oe  o'  N  oe  o  n 


-n 


nx 


=  Z 


n=l 


eNn  '  Cn 


(A-8) 
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and 


e0D= 


nx 

nx. 

nx 

E  0M 

n=l  n 

n=l 

nx 

E  .SNn  ' 

n=l 

Det 


(A-15) 


From  (A-9) ,  (A-10) ,  (A-14)  and  (A-15) ,  then  one  has 


1  f  nx  /  nx  \  nx 

8oe  =  S0  +  Det  z  5Nn  •(  £  Vftn  '  £  tn 

I  n=l  \n=l  J  n=l  I 

1  f  nx  .  /  nx  \"1 

**  ’  ^  '  1-^)1 


oe 


0o  + 


(A-16) 


(A- 17) 


(A-16)  and  (A-17)  are  results  for  any  specific  replication. 
Let  TF  be  the  time  from  first  measurement  to  fire,  then 


eeTF  =  *oe  +  9oe  •  TF 


Substitute  (A-16)  and  (A-17)  into  (A-18) ,  one  has 


(A-18) 


a  =0+6  *  TF  + 

8  eTF  o  o 


1  nx  l/nx  \ 

^  n=iN«P-lCV 


+  (tn  •  nx 


nx  ”| 

'  n-ltn>  •  TF| 


nx 

t  ill 

n  n=l 


(A-19) 


Define  E6  ,  the  error  ratio,  as  follows: 


11 


Efi  =  (fiP  -  e)/0  (A-20) 

Since 

6P  =  0 cTF  (A-21) 

and 

§  =  6  +  0  •  TF  (A-22) 

o  o 


Substitution  of  (A-21)  and  (A-22)  into  (A-20),  yields 


E0  =  0 


eTF 


-  ($0  +  eo  •  TF> 


9 


Also  equation  A-19  can  be  written  as 


(A-23) 


1  nx 

®eTF  '  ^ o  +  eo  '  TF)  “  Det  1  ® 

n=l 


(nx  »  nx  r  nxsT 

£  t„2|  -  t  •£  t_  +  f  t  • nx  -  £  t_\TF 

n-l  n)  n  n-1  "  \n  n-1  7  I 


Substitute  (A-24)  into  (A-23),  one  obtains 


nx 

E0=  Detn 


Knx  \  nx  ✓  nx  v  1 

s-M**  ’  nx~u  Sn 


(A-24) 


(A- 25) 


If  one  lets 


t 


<nx  \  /nx  \  /  nx  \ 

\wr  ' nx  ■  (A-26) 


and 

» 


* 


12 


(A-13) 


nx  /  nx  \ 

Det  =  nx  •  E  t  2  -(  E  tn ) 2 

n-1  n  V  n-1  7 

Substituting  (A-26)  into  (A-25) ,  yields 


E6  = 


nx 

e 

n=l 


^n 


6 


(A-27) 


Equations  (A-27),  (A-26),  and  (A-13)  are  similar  to 
aquations  (2),  (3),  and  (4)  in  the  text. 
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APPENDIX  B.  DERIVATION  OF  ERROR  IN  UNE- 
OF- SIGHT  RA3E  OF  ACCELERATION 

Similar  to  appendix  A  and  by  using  equation  (A- 17)  in  appendix  A 


■  nx 

1  Z  #  nx 

0=0  +  Det  .  n=l  §»,  (t  '  nx  -  Z  t  )  . 

oe  o  Nn  n  n 

Define  E0,  the  error  ratio,  as  follows: 

E0  =  (0P  -  0)/0 


also 


and 


0P  =  0 


0=0 


oe 


(A-17) 


(B-l) 


(B-2) 


(B-3) 


The  line-of-sight  acceleration  is  assumed  to  be  constant  with 
time.  Substitute  (B-2)  and  (B-3)  into  (B-l).  Therefore, 


E0  =  (0  -  0  )/0 

v  oe  o 


(B-4) 


and  equation  (A-  17)  can  also  be  written  as 


oe 


0  =  1 
Det 


nx  #  nx 

E  eK,  (t  *  nx  -  z  t  ) 
n.!"11  n-1  n 


(B-5) 


Substitute  (B-  into  (B-  4") ,  one  obtains 


E0  =  Det 


nx 


nx 


l  e^.  (t  nx  -  Z  t  ) 
n-l  141  n  n-1  n 


(B-6) 


Let 


nx 

CTDD  =  t  ‘  nx  -  Z  t 
n  n  n=l  n 


(B-7) 


Substitute  (B-7)  into  (B-6),  one  obtains 
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Ee 


_i_  ?x 

9 • Det  n=l 


•  CTDDn 


(B-8) 


Note  that  equations  (B-7)  and  (B-8)  are  similar  to  equations 
28  and  29  in  the  main  part  of  the  report. 


LIST  OF  SYMBOLS 

a  Angle  between  VT  and  line-of- sight  from 

launcher  (stationary)  through  the  target 
at  time  of  fire  or  X  axis. 


a 


t 


6 

• 

0 

0 

• 

0 

0 


en 

eTF 

Mn 

Nn 


0 

0 


0 


oe 


0 


0 

o 

®  oe 
A 

AP 

C  and  CTDD 
n  n 

Det 


Angle  between  acceleration  direction  of 
target,  at  time  of  fire,  and  X  axis. 

Line  of  sight  rate  at  time  of  fire. 

n-th  value  of  0,  estimated. 

0  estimated  at  time  of  fire,  TF. 

n-th  value  of  0,  measured. 

0  noise  contribution  for  the  n-th 
measurement . 

True  initial  value  of  0  at  time  of  first 
measurement . 

Initial  estimated  value  of  0  of  target. 

Line-of- sight  acceleration  rate  (first 
derivative  with  time  of  §). 

n-th  estimate  of  0 . 

True  value  of  0  of  target. 

Estimated  value  of  0  of  target. 

Acceleration  of  target. 

Angle  between  velocity  of  projectile  and 
X  axis. 

Defined  in  equations  3  and  29. 

Value  of  determinent,  defined  in  equation 
4. 


Dt  Time  between  each  measurement. 

E0  Error  ratio  in  0  (equations  2  and  8) . 
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Ee 

Error  in  e  double  dot,  divided  by  0 
(equations  9  and  28). 

I 

Denotes  component  in  direction  of  X 
axis . 

J 

Denotes  component  perpendicular  to  X 
axis . 

Ki,K2 

Constants  defined  in  equations  25  and 

26. 

MD 

Miss  Distance  between  projectile  and 
launcher  (in  direction  perpendicular 
to  X  axis) . 

fMAX  or 
nx 

Maximum  number  of  measurements . 

P 

Perceived  value  of  the  term  preceding 

P;  perceived  by  the  launcher  either 
from  measurement  or  the  launcher's 
calculations . 

SIGMA  MD 

Sigma  miss  distance. 

t 

n 

Time  of  the  nth  measurement . 

TF 

Time  from  first  measurement  in 
prelaunch  phase  to  fire. 

TINT 

Time  to  intercept. 

VP 

Velocity  of  projectile. 

VT 

Velocity  of  target. 

VTJ 

Velocity  of  target  in  the  J  (perpendicular) 
to  X  axis)  direction. 

X 

0 

Distance  from  launcher  to  target  at 
time  of  fire. 
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